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Abstract 

Using recent results for nonperturbative contributions to the B and D meson inclusive 
semileptonic widths, a model independent extraction of \Vhc\-, rric and mb is made from the 
experimentally measured B and D lifetimes and semileptonic branching ratios. Constrain- 
ing the parameters of the HQET at 0{l/mQ) by the D semileptonic width, \Vi,c\ is found 
to lie in the range .040 < iV^d < 0.057. The c and b quark masses are not well constrained 
due to uncertainty in the relevant scale of as- These results assume the validity of pertur- 
bative QCD at the low scales relevant to semileptonic charm decay. Without making this 
assumption, somewhat less stringent bounds on Vbc from B decay alone may be obtained. 



UCSD/PTH 93-25 
UTPT 93-21 
CMU-HEP 93-13 

DOE-ER/40682-37 Revised version, October 1993 



1 



There has been much recent interest in the apphcation of the techniques of the heavy 
quark effective theory (HQET) to inclusive decays of hadrons containing a single c 

or b quark [H-fl^. As first made explicit in [0, the differential decay rate for an inclusive 
factorisable process such as B ^ XcCVe or B ^ Xg'j may written, using an operator 
product expansion, in terms of the matrix elements of local operators between B mesons. 
Using the techniques of HQET, it was shown in that the leading term of the expansion 
in l/mh reproduces the parton model, and that the 0{l/mi,) corrections vanish by the 
equations of motion for the heavy quark. More recently, the 0{l/ml) corrections to the 
parton model have been calculated for semileptonic decays fl^lQ? as well as for the rare 
decays B Xg'j and B Xse^e~ ||12|| . A crucial observation which emerges 

from the operator product expansion is that the semileptonic width of a heavy meson is 
proportional to the fifth power of the heavy quark mass, not the meson mass. The heavy 
quark mass nig is a well-defined quantity in HQET, and is defined such that the residual 
mass term in the effective heavy quark Lagrangian vanishes (for a detailed discussion of 
the definition of ruq, see [1T3[|). 

In addition to these nonperturbative corrections, there are calculable (9(as(mg)) cor- 



rections to the free quark decay picture coming from real and virtual gluon emission |n 
In this letter, the complete expression to O (as(mg), for the semileptonic width 

of a heavy meson is used to extract the c quark mass from the measured lifetime 
and semileptonic branching ratio. From this, the b quark mass is determined and the 
B lifetime predicted in terms of the KM matrix element \Vbc\- Comparing this with the 
experimentally measured lifetime allows |Vbc| to be determined. 

Throughout this letter the mass of a heavy meson will be denoted by Mq and the 
mass of the corresponding quark by m^. The meson and quark masses are related by [^| 



The parameter A in (P may be interpreted as the energy of the light degrees of freedom of 
the meson. It corresponds to the constituent mass of the light quark in the meson, and so is 
expected to be of order a few hundred MeV. From consideration of QCD mass inequalities 
analogous to those used in the light quark sector to prove that nip > m^^, Guralnik and 



Manohar IjT^] have shown the existence of a rigorous lower bound 

A > 237 GeV. 
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We will also take A < 800 MeV in this work when quoting limits on nic, rrii, and Vbc- This 
bound is consistent with the QCD sum rules estimate (see also the second reference 
mil) 

a'''' = 570 ±70 MeV (3) 

and is also consistent with the quark model, in which A is expected to be a few hundred 
MeV. 

Ai and A2 are defined in terms of the expectation values in the heavy quark effective 
theory, 

< HQ\h{iDfh\HQ >=2MqXi 

(4) 



< HQ\h a'^'G^^ h\HQ >= 6MqA2(^) 

where Hq is the pseudoscalar heavy meson {B or D) and h is the heavy quark field in 
the effective theory. A2 parameterises the effects of the chromomagnetic moment operator, 
and may be extracted from the B* — B mass splitting 



A2(m;,) = ^(Ms* - Mb) ^ -AMI, - M|) = 0.12 GeV^ 
o 4 



(5) 



(corresponding to A2(toc) = 0.10 GeV^ HI])- Ai parameterises the kinetic energy of the 
b quark inside the hadron. Since it contributes equally to and m^*, Ai cannot be 
extracted from the meson masses. 

At (9(l/mg), the nonperturbative corrections to the parton model are also parame- 
terised in terms of Ai and A2. Combining the results from [0-|ll[] with the as corrections 
to the free quark model [113 gives the complete expression for the semileptonic D width to 
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The functions /i,2 arise from the finite mass of the quark in the final state, 

/i (x) = 1 - Sx^ + 8x^ - x^ - 24x^ log X 

8 5 
/2(x) = 1 - -x^ - 8x* + 8x^ + -x^ + 8x^logx, 
3 3 



(7) 



while the function g{x) arising from one-gluon graphs is tabulated in Ref. (we note 
thatt/(0) = 3.62, ^(.2/1.6) = 3.15 and ^(1.6/4.8) = 2.40). We also take = 200±80MeV 
and nid — 0. 

The choice of the appropriate scale fXc in is a potentially troublesome one. The 
scale is set roughly by rric; however, since the final state hadron typically carries off only 
a fraction of the available energy, the appropriate scale may be somewhat lower. Formally 
this choice of scale is a higher order effect; however, in practice there is a considerable 
difference between, for example, as{rnc) and as{mc/S). In particular, at scales much less 
than mc, QCD perturbation theory cannot be trusted. Since we have not done the two-loop 
calculation, we cannot resolve this scale ambiguity. In this letter we will take fic — i^c, 
and the sensitivity of the results to higher order terms in ag will be estimated by formally 
varying /ic between rric and rric/S. Taking Aq^q = 250 MeV and Aq^q = 300 MeV, and 
using the two-loop expression for as{n) gives as{mc) = 0.09 tt and as{mc/^) — 0.4 tt. We 
stress that we are not claiming that QCD perturbation theory is valid at /U ~ rric/ 3, but 
rather that if perturbation theory is meaningful for D decays, we expect this approach 
to give a reasonable estimate of the uncertainty due to higher order QCD corrections.^ 
Fortunately, the determination of Vbc is relatively insensitive to the uncertainty in rric 
extracted from the D width. 

From the observed semileptonic branching ratio, lifetime and mass of the ETH 

Br(D± ^ + anything) = (17.2 ± 1.9)% 

T^± = 10.66 ±0.23 X 10"^^ s (8) 
Mo± = 1869.3 ±0.3 MeV 

it is a simple matter to extract rric and Ai from Eqs. (^ and (P for any value of A. (We 
have checked that the gives results consistent with the D^; the meson is used, as 
the uncertainty in its semileptonic branching fraction is somewhat smaller than for the D^.) 
The results are shown in fig. |l| (a) and (b). In each figure the narrow band corresponds 
to He = rric, while the broader band correspond to rric/S < Hc < rric- The uncertainty 
corresponding to la errors on all experimental input parameters is included in the width 
of each band. From fig. |1|, the Guralnik-Manohar bound (^ on A translates to limit 

Ai > -0.5GeV (la). (9) 



^ Note that the energies encountered in this problem are similiar to those encountered in 
hadronic r decays, for which the corrections to O (as(?T^T)^) have been calculated, and for which 
the perturbation series appears to be converging ||20||. 
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Because of the uncertainty in |Uc, the c quark mass is not well constrained from the data, 
although taking A < 800 MeV suggests rric > 1460 MeV. Since the upper limit of the error 
bar in fig. [^(a) is both very sensitive to the value of as at low scales and corresponds to 
rric > mo, a useful upper bound on rric cannot be extracted. 

For B mesons, formulas analogous to (|l|) and hold, with the replacements Mo — ^ 
Mb, rric T^h, Vcs Vbc, Vcd Vbu and m^. — > rric. Since Ai is now determined from D 
decays for each value of A, the b quark mass may be solved for as a function only of 
A. Using Mb = 5.278 ± 0.002 GeV gives the results shown in fig. 0(c). Since Mb is less 
sensitive to Ai than Md, it is not surprising that rrib is roughly a linear function of A, with 
smaller error bars than rric. i^b is found to lie in the range 5140 MeV > > 4600 MeV. 

Given rrib and Ai, the semileptonic width of the B may now be predicted, allowing 
the KM matrix element \Vbc\ to be determined as a function of A. The inclusive branching 



ratio of the B to electrons is measured to be 21 



Br(S ^ e + anything) = 10.7 ± 0.5% (10) 
(where the charge of the B is not determined). Averaging the recent measurements of the 



B lifetime by the DELPHI collaboration and by the CDF collaboration gives a B 
lifetime of 

tb = 1.32 ±0.09 X 10"^^ s (11) 

(statistical and systematic errors added in quadrature). This gives \Vbc\ as a function 
of A, as shown in fig. ^ The two bands in the figure correspond to the two choices of 
scale fXb = rnb, lie = rUc and rric/S < fXc < rn^ rrib/S < /ib < rrib (with fXc/ fJ-b fixed). 
Over the range of values for A shown in fig. |^, \Vbc\ is found to lie in the range 0.040 < 
\Vbc\ < 0.057. The lower bound is set by the rigorous lower bound on A, while the upper 
corresponds to A = 800 MeV. This extraction is consistent with the current value of 

1/2 

\Vbc\ (1.32 ps/rs) ' = 0.051 ±0.008 extracted from the exclusive decay B D*eVe using 
HQET Jl^ . We note that choosing a lower scale for Hc and Hb tends to raise the preferred 
values of \Vbc\. 

For a given value of A the effects of higher order terms in l/rriq and ag on this 
extraction of \Vbc\ may be estimated. These terms may be parametrized by db,c and e^^c, 
defined by 

A ^ Ai±3A2 , 
Md =mc 1 ± ' ^ + 5, 
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They correspond to a change in \Vbc\ of 



\Vbc\ 




Set + 5(5, -35,)). 



(14) 



The leading contributions to Sc are of order (A/rric)^ and ""^J^"-^ A2/?Tie^, which we estimate 
to be at most a few percent, while 5, is expected to be at least an order of magnitude 
smaller. We have attempted to take the leading contributions to €,,6 arising from higher 
order QCD corrections into account by varying the scales fj,c,b between mb,c and mb,c/'i- 
From dU), we can understand why our results were relatively insensitive to the choice 
of scale: an uncertainty in the charm width only gives ~ 1/6 this uncertainty in \Vbc\- 
Furthermore, the uncertainty in scale tends to cancel between and e,. 

Since the extraction of Ai as a function of A depends on the validity of perturbation 
theory for QCD at a low < rric) scale, it is instructive to compare this extraction of 
\Vbc\ with that obtained by ignoring the D decay data altogether and simply varying Ai 
between ±1 GeV^ for each given value of A. This amounts to working at O (1/mg) in the 
heavy quark expansion with a reliable estimate of the uncertainty from l/m^ terms. As 
before, the scale fib is also varied between fXb = rub/^ and fib = rnb- The results are shown 
in fig. 1^. The shaded area on each graph corresponds to both the experimental errors on 
the input parameters and the variations in Ai and fib- As expected, the uncertainty is 
\Vbc\ is somewhat larger than when the D decay data are included, particularly at large 
values of A; however, \Vbc\ is still well constrained since the prediction for the B width is 
relatively insensitive to Ai. For 237 MeV < A < 800 MeV, \Vbc\ is found to lie in the range 



It is encouraging that the extraction of \Vbc\ from inclusive decays is consistent with 
that obtained from exclusive decays over the full allowed range of A. Furthermore, the 
extraction of \Vbc\ from inclusive decays has some theoretical and experimental advantages 
over exclusive decays. No exclusive hadronic final state needs to be identified, and the 
extrapolation of form factors to zero recoil is not required. The method is limited by 

^ The operator h(iD)^h receives no strong interaction corrections to its coefficient, by repa- 
rameterisation invariance [p4| . 



0.04 < \Vbc\ < 0.06. 
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the uncertainties in A and the relevant hadronic scales fit, and fic- Additional input, 
such as a lattice measurement of A or Ai, or a two-loop calculation of the semileptonic 
decay rate, would further constrain \Vbc\- For example, from fig. |l](a), the QCD sum 
rules estimate (D of A = 570 ± 70MeV corresponds to 0.1 GeV^ < Ai < 1.5 GeV^ and 
0.044 < |Vf,c| < 0.054 (although the upper bound on Ai should not be taken too seriously, 
since both the l/rUc expansion and QCD perturbation theory are unreliable at this point). 
On the other hand, the nonrelativistic quark model suggests that A ~ 350 MeV, the mass 
of the light constituent quark, and that Ai ~ — (350MeV)^ ~ — O.lGeV^ (note that Ai < 
corresponds to a positive mass shift from the kinetic energy of the heavy quark). From 
fig. 0(a), these values of Ai and A are consistent with one another, and are consistent with 
somewhat lower values of \Vbc\- 

We thank M. B. Wise for useful discussions, particularly concerning the relevant scales 
fXc and and the possibilities of extracting Vbc without using the D decay data. We are 
also grateful to A. Falk, B. Holdom, A. Manohar and L. Wolfenstein for useful comments 
and criticisms. MJS acknowledges the support of a Superconducting Supercollider Na- 
tional Fellowship from the Texas National Research Laboratory Commission under grant 
FCFY9219. This research was supported in part by TNRLC grant RGFY93-206 and by 
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Figure Captions 



Fig. 1. Ai, rric and as functions of A. In each graph, the vertical hne corresponds to 
the Guralnik-Manohar bound A > 237 GeV. The shaded regions correspond to 
the experimental error in the input parameters and, for the broader band, the 
variation in Hc- 

Fig. 2. The weak mixing angle iV^d as a function of A for tb = 1-32 ± 0.09 ps. The 
vertical line corresponds to the constraint imposed by the bound on A. The error 
bar on the left indicates the current value of \Vbc\ extracted from the exclusive 
decay B D*eVe and its associated uncertainty. The shaded region corresponds 
to the experimental error in the input parameters and, for the broader band, the 
variations in Hc and Hb- 

Fig. 3. The weak mixing angle \Vi,c\ as a function of A for tb = 1-32 ± 0.09 ps, without 
imposing constraints on Ai from D decay. The vertical line corresponds to the 
lower bound on A, and the shaded regions correspond to the experimental error 
in the input parameters, varying Ai from — IGeV to 1 GeV and fii, from nii, to 
mb/3. The error bar on the left indicates the current value of |V5c| extracted from 
the exclusive decay B — > D*eVg and its associated uncertainty. 
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Figure 1 (b) 
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Figure 1 (c) 
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